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Abstract — This paper considers the input-erasure Gaussian 
channel. In contrast to the output-erasure channel where 
erasures are applied to the output of a linear time-invariant 
(LTI) system, here erasures, known to the receiver, are applied 
to the inputs of the LTI system. Focusing on the case where 
the input symbols are independent and identically distributed 
(i.i.d)., it is shown that the two channels (input- and output- 
erasure) are equivalent. Furthermore, assuming that the LTI 
system consists of a two-tap finite impulse response (FIR) 
filter, and using simple properties of tri-diagonal matrices, an 
achievable rate expression is presented in the form of an infinite 
sum. The results are then used to study the benefits of joint 
multicell processing (MCP) over single-cell processing (SCP) in 
a simple linear cellular uplink, where each mobile terminal is 
received by only the two nearby base-stations (BSs). Specifically, 
the analysis accounts for ergodic shadowing that simultaneously 
blocks the mobile terminal (MT) signal from being received 
by the two BS. It is shown that the resulting ergodic per-cell 
capacity with optimal MCP is equivalent to that of the two- 
tap input-erasure channel. Finally, the same cellular uplink is 
addressed by accounting for dynamic user activity, which is 
modelled by assuming that each MT is randomly selected to be 
active or to remain silent throughout the whole transmission 
block. For this alternative model, a similar equivalence results 
to the input-erasure channel are reported. 

I. Introduction 

In recent years various variants of erasure networks have 
been the focal point of many works in the fields of infor- 
mation theory, coding theory, communication, and digital 
signal processing (DSP). It turns out that the idealization 
of situations where symbols are detectable with very high 
or very low probability is useful in providing accurate 
and tractable assessment of practical interest. Applications 
of erasure networks, in general, and erasure channels, in 
particular, include packet-switched networks, wireless fading 
channels, impulse noise channels, and sensor networks with 
faulty transducers. Another interesting DSP application is 
decimation using frequency shaped random erasures [1] (see 
also [2] for a somewhat dual problem). 

The simplest non-trivial setting incorporating erasures is 
perhaps a memoryless discrete channel observed through a 
memoryless erasure channel. In this case, it is straightforward 
to see that the capacity of the concatenated channel is 
the capacity of the discrete channel times the erasure rate 
[3]. Moreover, even if the erasures have memory, the latter 
result holds [3] [4]. However, adding memory to the discrete 
channel changes the picture considerably. 
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Fig. 1. The linear input- and output-erasure Gaussian cliannels. 



In [5], Tulino et al. have calculated the capacity of an input 
power-constrained Gaussian channel concatenated with a 
memoryless erasure channel. The problem requires obtaining 
the asymptotic spectral distribution of a sub-matrix of a 
nonnegative definite Toeplitz matrix obtained by randomly 
deleting rows and columns. In addition, [5] presents an 
explicit expression for the optimal input power spectral 
density (PSD) function. Asymptotic capacity expressions in 
various regimes such as extreme-SNR, sporadic erasures 
events, and sporadic non-erasure events, are reported as well. 
In addition, several useful bounds (lower and upper) in closed 
form are presented. Recently, the results of [5] were extended 
in [6] to include any flat fading distribution, not necessarily 
the on-off process obtained by erasures. 

In this paper, we consider the input-erasure Gaussian 
channel. According to this setup, erasures (known at the 
receiver only) are applied directly to the Gaussian input 
of a linear time-invariante (LTI) system (e.g. a finite input 
response filter, FIR), whose output samples are affected 
by additive Gaussian noise. To simplify the analysis, we 
consider independent and identically distributed (i.i.d.) in- 
put signaling, and show that in this case the input-erasure 
channel and the output-erasure channel of [5] are equivalent. 
Focusing on the case where the channel filter includes two 
taps, we present an alternative analysis to [5] that obtains 
the achievable rate without relying on random matrix theory 



tools. The analysis is based on a simple observation that in 
the two-tap case, each erasure splits the output vector into 
two independent sub-vectors. In addition, it uses a useful 
property of tri-diagonal matrices, which states that their 
determinant can be expressed as a weighted sum of the 
determinants of their two main principal sub-matrices. It is 
also shown that the analysis can be extended to the case in 
which the erasures are not memoryless. 

The results are then used to find the achievable rate 
supported by a simple cellular uplink model (where inter- 
cell interference is limited only to adjacent cells), in which 
the transmission of each mobile terminal (MT) undergoes 
erasures due to ergodic shadowing, which simultaneously 
blocks it from being received by all the base-stations (BSs). 
We show that with optimal joint decoding of the received 
signals from all the BSs (that is, multicell processing, MCP), 
the ergodic per-cell capacity of this cellular uplink equals the 
rate of the two-tap input-erasure channel considered earlier 
The performance of MCP is compared to that of single-cell 
processing (SCP), demonstrating that the benefits of MCP 
decrease with the erasure rate (the reader is referred to [7] 
and [8] for recent surveys on MCP). We then consider a 
similar cellular uplink setup but with a dynamic user activity 
model, in which each user is randomly active or silent for 
the whole transmission block. We show that with MCP the 
system throughput per active user is equal to the rate of 
the two-tap input-erasure channel divided by the probability 
of an arbitrary user being active. Similar benefits of MCP 
over SCP, as those reported for the shadowing scenario, are 
observed also for the user activity model. 

It is finally noted that a related analysis has been recently 
reported in [9], where the outage capacity supported by a 
similar cellular uplink and user activity model was derived 
with both MCP and SCP 

II. Channel Model 

In this work we consider a channel (see Fig. [T]a), where 
erasures {e„}, known to the receiver only, are applied to 
the inputs symbols of a discrete time LTI system with 
transfer function H{f). For simplicity, we assume that the 
LTI system is causal and has a finite^ impulse response of 
length [L + 1) symbols. Accordingly, the received signal at 
time index n, is given by 



m=0 



(1) 



where {^m}m=o '■^^ impulse response coefficients such 
that 

L 

HU) = E , (2) 

m— 

and Zn is the additive complex Gaussian noise z„ ~ 
CA/'(0, 1). The erasures e„, are assumed to be i.i.d. Bernoulli 

'it is easily verified tliat tlie results extend to any LTI system with finite 
g?:mJ^\H{f)\'df<cx,. 



distributed random variables (rv.'s) e„ ~ B{q) (i.e., e„ S 
{0, 1}, Pr(e„ - 0) = g, Fr(e„ = 1) = 1 - q). 

We consider transmissions in blocks of N symbols. The 
received vector of an arbitrary block is given by 



Y, 



N+L - 

where Yn+l = (yi, 



(3) 



jUn+lV is the {N + L) x 1 output 
vector, Xn — (a;i, . . . , xjv)^ is the complex Gaussian x 1 
input vector ^ CJ\f{0, Qjy) with average power con- 
straint tr (Qjv) < XP, En = diag(ei, . . . , cn) is the N x 
TV diagonal erasure matrix, and Z]y-^-L — {zi, ■ ■ ■ , zat+i)^ is 
the (A^ + i) X 1 noise vector Zn+l ~ CAA(0, Jjv+l)- In 
addition, Hn+l is the {N + L) x N L-diagonal Toeplitz 
channel transfer matrix with [HN+L]n,m — hn-m (where 
out-of -range indices should be ignored). 

Recalling that the receiver is aware of the erasures, it is 
easy to prove that a Gaussian input distribution is a capacity- 
achieving distribution and that the capacity is given by 



C 
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tr(Q^)<NP 
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(4) 



where the expectation is taken over the diagonal entries of 
the erasure matrix En, and In+l is an (N + L) x {X + L) 
identity matrix. 

Unfortunately, finding the optimal covariance matrix and 
calculating (HJi is still an open problem. Instead, we focus in 
the rest of this work on assessing achievable rates resulting 
by using i.i.d. inputs, hence, by setting ~ PIn- 

III. Rate Analysis 

A. Arbitrary FIR Filter 

Here we consider an arbitrary L-length FIR filter and 
rewrite the logdet term of ^ as: 



logdet (In+l + PHn+lEnE^^H' 



log det ( / 
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N 



(5) 



By noting that the left-hand-side (LHS) of the last equation 
is the same as the one resulting from the analysis of the 
output-erasure channel with i.i.d. inputs, considered in [5] 
(see Fig. [T]b), we can state the following. 

Proposition 1 The Gaussian input-erasure and output- 
erasure channels with i.i.d. inputs are equivalent. 

This result implies that with i.i.d. inputs, both channels 
achieve the same rate. Hence, all the results reported in [5] 
hold verbatim for the input-erasure channel with i.i.d. inputs, 
by setting the PSD accordingly (i.e., S{f ) = P \H{f )\^) in 
the results of [5]. 

The analysis in [5] uses random matrix theory tools to 
derive the main result which is the rate of the output-erasure 
channel for any channel FIR filter and any input PSD. The 
rate result of [5] involves a fixed-point integral equation and 



in general it is not explicitly formulated. Next, we present 
an alternative analysis for the special case of i.i.d. inputs 
and a two-tap channel filter. The resulting rate expression is 
considerably different than that of [5], and may shed some 
further light on the problem. 

B. Two-Tap Filter (L — 1) 

In this section we consider a two-tap channel filter (i.e., its 
impulse-response includes only two non-zero coefficient ho, 
and hi). Hence, the channel transfer matrix H m+i reduces 
to a bi-diagonal Toeplitz matrix. For this special case we 
have the following. 

Proposition 2 The achievable rate of the two-tap input- 
erasure Gaussian channel with i.i.d. inputs is given by 



i?2tap = g' 5^(1-9)" log 



r — s 



(6) 



where 



a^l + P{\hoY + \hi\') and b^P\ho\\hi\ , (7) 
and 

r = - + \J a?- — 46^^ and ~ 2 ("^ ^ \/ a? — 46^ 

(8) 

Proof: See Appendix ■ 
The proof relies on a simple observation that in the two- 
tap case, each erasure splits the output vector into two 
independent sub-vectors. In addition, it uses a useful property 
of tri-diagonal matrices, which states that their determinant 
can be expressed as a weighted sum of their first two main 
principal sub-matrices. It is noted that this analysis does not 
involve classical random matrix theory. Moreover, since it 
relies heavily on the two-tap assumption, it cannot be applied 
to input-erasure channels with filter lengths larger than two. 

Examining the rate expression (|6]l, a few comments are 
in place. Although the rate expression is an infinite sum, 
any finite sum, which is evidently a lower bound, can 
be calculated numerically in a straightforward manner. In 
addition, convergence is assured since the expression is a 
power series and (1 — g) < 1. It is also observed that 
the rate is independent of the filter coefficients' phases, and 
hence, random phases (known to the receiver) can be added 
to the channel model without changing the end result. By 
applying Hadamard's inequality directly to the covariance 
matrix {I m + Gn) (see (|27]i), it can be shown that for a 
fixed filter gain, filter memory reduces the rate. Hence, the 
rate is maximized for a one-tap filter with equal gain. On this 
note, it is easily verified that for the special case of /ii = 
(memoryless Gaussian channel with erasures), s = and (|6]l 
reduces to 



i?ltap = (l-g) log(l+P|;io 



(9) 



which was already reported in [10]. Another special case, 
which is also an upper bound, is the case where <? = 0. The 



capacity of this two-tap channel with i.i.d. inputs is given by 
[11] 

^2tap-ub = logr , (10) 

where r is defined in dHJ. 

Turning to asymptotic analyses, we refer the reader to 
[5], where compact closed-form expressions are reported for 
the rate low-SNR characterization, the rate with sporadic 
erasures, and the rate with sporadic non-erasures. These 
expressions can be used verbatim for the two-tap input- 
erasure channel of interest by setting i.i.d. inputs in [5]. 
In addition, [5] provides also the high-SNR characterization 
of the output-erasure channel rate, albeit the calculation of 
the high-SNR power offset involves a fixed point-equation. 
Applying definitions of [12] to the rate expression (|6]l, we 
provide an alternative expression for the high-SNR power 
offset as well. 

Proposition 3 The two-tap input-erasure Gaussian channel 
rate with i.i.d. inputs is characterized in the high-SNR regime 
by 



5*00 = (1 - q) and 

00 
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(11) 



And the deleterious effects of increasing erasure rate q are 
clearly visible. 

1 ) Erasures with Memory: Following the proof of Prop. 
121 it can be verified that the presented analysis method can 
handle also erasures with memory (stationary and ergodic). 
The only stage that should be altered is the calculation of 
dSST l, that is, the probability of having an isolated sequence 
i of n consecutive non-erasures should be recalculated for 
the new erasure source. For example, assuming the erasures 
follow a first order Markov chain with transition probabilities 
Pr{0 1) = 1 — qo, and Pr{l — > 0) = qi, it is easily 
verified that 



Pr{£- 



qK^-qo) 



a-qi) 



(12) 



1 - go + 91 

and the corresponding rate of the two-tap input erasure 
channel with i.i.d. inputs is given by 



J^2tap — 



g?(i-go) 

1- qo + qi 



fl-gi)"-Mog 



n=l 



„n+l 



r — s 

(13) 

As expected, for qo = qi = q the erasure source is i.i.d. 
and ( fTSl ) reduces to It is interesting to note that for the 
special case hi — 0, which corresponds to a memoryless 
channel followed by an erasure channel with memory, the 
rate reduces to 



i?ltap = (l-g)log(l + P|/lor) , 



where 



qi 



1 - go + 91 



(14) 



(15) 



Fig. 2. Rates vs. erasure rate q, for P = 0, 2, 4, 6, 8, 10, 12, and 
\hof = 0.8, = 0.2. 
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Fig. 3. Rates vs. SNR P, for q = 0, 0.05, 0.1, 0.2, 0.3, 0.4, and 
|/io|^ =0.8, = 0.2. 



is the steady-state erasure rate. This results is in agreement 
with [4], where it has been shown that the capacity of a 
concatenation of a discrete memoryless channel with capac- 
ity C, and an erasure channel (possibly with memory) with 
erasure rate q, is equal to {1 — q)C. 

IV. Numerical Results 

The rate ^ is plotted in Figure |2] as a function of the 
erasure rate q, for several SNR values, and filter coefficients 
|ft-ol^ = 0-8, 1^1 1^ = 0.2. Expression (|6]l is calculated 
by taking 200 summands and is plotted with solid lines, 
while Monte-Carlo simulation results of 50 trials, each with 
block size of 200 symbols, are shown as well. Examining 
the figure, the deleterious effect of erasures is evident. In 
addition, the figure demonstrates a good match between the 
Monte-Carlo simulations and the "exact" results. Similar 
conclusions can be obtained from Figure [3] which includes 
rate curves (calculated from ^ as above) as functions of 
the SNR P for several values of the erasure rate q, and 



Fig. 4. Rates vs. filter coefficient j/ioP. for 9 

0, 0.05, 0.1, 0.2, 0.3, 0.4. 



ft-ol^ — 0.8, = 0.2. In addition, the erasure-free 

(q = 0) rate curve ( fTOl i is included as a reference. 

In Figure |4] the rate curves are plotted as functions 
of the filter coefficients \ho\^ for a unit gain filter (i.e., 
\ho\^ + \hi\^ = 1), and P = 10 [dB]. It is observed 
that the worst-case filtering corresponds to coefficients with 
equal amplitudes (i.e., |/io|^ = |/ii|^). On the other hand, 
as claimed earlier, best-case filtering is achieved for flat 
or memoryless channels (i.e., \hi\^ = or \hi\^ = 0). 
Moreover, these observations are independent of the erasure 
rate q. 

V. Application to Cellular Communications 

Cellular systems are currently the main media to pro- 
vide high-data rate services to mobile users. Therefore, a 
continued search for techniques that provide better service 
and coverage in cellular systems is under way. A promis- 
ing cooperation-based technique is joint MCP, proposed by 
Wyner in [13]. Accordingly, clusters of BSs jointly process 
their signals to mitigate or eliminate the overall interfer- 
ence, since more signals are useful. Unlike conventional 
approaches that treats interference stemming from other cells 
as noise, or tries to avoids it, this approach (also referred 
to as distributed antenna array) exploits the interference to 
provide better rates to the clusters' MTs. 

Here we are interested in further exploring the benefits of 
MCP when ergodic erasures affect the MTs' transmissions 
in the uplink. Such erasures account for a scenario where an 
MT may undergo ergodic deep fading (or shadowing) such 
that its transmission cannot be received by any of its cluster's 
antennas. 

Motivated by the fact that inter-cell interference is es- 
sentially limited to only a small number of BSs, and with 
mathematical tractability in mind, we assume a synchronous 
Wyner-like cellular uplink setup. According to this setup 
(see Fig. |5]l, an infinite number of cells are arranged on a 
line, as along a highway. Moreover, the each MT "sees" 
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Fig. 5. The soft-handoff cellular uplink model. 



its own BS antenna (with unit path gain), and the adjacent 
BS antenna (with path gain a G [0, 1]) only. This model, 
which was first introduced in [11], focuses on users lying 
on the cells' boundaries and is thus sometimes referred to 
as the soft-handoff model. Specifically, assuming a single 
active user per-cell (intra-cell time-division multiple-access), 
the received signal at the nth BS antenna, for an arbitrary 
time index, is given by 



Vn 



(16) 



where Xn and are the MTs' transmissions Xn, Xn-i 
CM{0,P), e„ and e„_i are the erasures e„, e„_i ^ B{q), 
and Zn is the additive noise z„ ^ CAf{0, 1). It is noted that 
the erasures are assumed to be i.i.d. among different users 
(along the cells' indices n), and ergodic along the time index 
for each user. Finally, it is noted that users are not allowed 
to cooperate. 

A. Single Cell Processing 

As a reference results we consider two conventional SCP 
approaches. 

1) No Frequency Planning: The simplest approach is for 
each cell to decode its own signals while treating the signals 
from the interfering cell as noise. Assuming that each cell is 
aware of its own user's erasures and those of the interfering 
cell's user (see [14] for a similar setting), it is easily verified 
that the achievable rate is given by 



-Rs. 



E 



log 1 



1 + a'^CintP 



(17) 



where the expectation is taken over Cgig and eint which 
represent arbitrary i.i.d. erasure r.v.'s of the useful signal and 
interference, respectively, egig, Cint ^ By straightfor- 

ward calculation we get that ([TT] ) is given by 



Rscp = (1 - g)log 



(l + F)'?(l + (l + a^)P; 



1-9- 



The high-SNR characterization of i?scp is given by 
Soo = q(1 - q) and Coo = . 



(18) 



(19) 



It is interesting to observe that the rate above demonstrates 
a non-interference limited behavior, although the receiver 
treats the other cell's signals as noise. This is easily explained 
by noting that with probability q{\ — q) the useful signal is 
received without interference. 

2) Inter-Cell Frequency Sharing: Here we consider a 
more involved approach that requires static resource shar- 
ing between the cells. Since, according to the soft-handoff 
model, interference stems from one cell only, by dividing 
the available bandwidth of W [Hz] into two orthogonal 
subbands with bandwidth of W/2 [Hz] each, and assigning 
them alternately to even and odd indexed cells, interference 
is totally avoided. Accordingly, it is easily verified that the 
achievable rate is given by 

1„.. l-q- 



i[log(l + eP)] 



log(l + P) , (20) 



where e is an arbitrary erasure rv., e ^ B{q). 

The high-SNR characterization of i?icfs is given by 

1-q 
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(21) 



where the deleterious effects of both the erasures, and the 
bandwidth sharing are clearly observed in the multiplexing 
gain expression. Comparing the multiplexing gains of the 
SCP and inter-cell frequency sharing (ICFS) rates, reveals 
the superiority of the latter for erasure rates g < 1/2 in the 
high-SNR regime. 

B. Multicell Processing 

With MCP, the BSs send their received signals to a central 
processor (CP) via an ideal backhaul network. The CP 
collects the received signals and jointly decodes the MTs' 
messages. Since the CP is aware of all the erasures, the 
overall channel is an ergodic multiple-access channel (MAC) 
for which the sum-rate capacity divided by the number of 
cells is given by 



Cmcp = lim 



E 



^ logdet [In+1 + Pi?N+l£^Ar-E]vi?jv+i 



(22) 



where En = diag(ei, . . . , bat) is the N x N diagonal 
erasure matrix, and Hn+i is an {N + 1) x bi-diagonal 
Toeplitz channel transfer matrix with [Jfjv+i]n.n — 1 and 
[HN+i]n,n-i = OL (where out-of-range indices should be 
ignored). It is noted that, by the symmetry of the model, the 
ergodic per-cell sum-capacity derived above is also the equal 
per-cell capacity. 

Comparing the rate expression of the two-tap input-erasure 
channel with i.i.d. inputs (see (|4|l) to the per-cell sum- 
rate supported by the soft-handoff model, the following is 
evident. 

Corollary 4 The per-cell sum-rate supported by the soft- 
handoff model with optimal joint decoding and no MT 
cooperation equals the rate of the two-tap input-erasure 
Gaussian channel with i.i.d. inputs and /iq = 1, hi = a. 



Fig. 6. Rates vs. erasure rate, for = 0.5 and P = 14 [dB]. 



Fig. 7. Throughputs per active user vs. non-active probability q, for oP = 
0.5 and P = 14 [dB]. 



Since the shoft-handoff model is equivalent in terms of its 
per-cell sum-rate to the input-erasure channel, all the results 
and conclusions reported in Section |III] also hold here. In 
particular the high-SNR characterization reduces to 

5*00 = (1 - q) and 

r yn-u (23) 

n— 1 ^ ^ 

The benefits of MCP are evident in the high-SNR regime 
since its rate has two-fold degrees of freedom relative to the 
ICFS rate, and l/q-fold degrees of freedom relative to the 
SCP rate. 

It is noted that the MCP rate is an increasing function of 
the inter-cell interference factor a. This is since: (a) r and s 
are increasing with a ; (b) r > s ; and (c) by recalling that 

n+l _ n+1 " 
^ y j.n-m^n, ^ ^24) 

r — s ^-^ 

m=0 

which is a summation of multiplications of increasing func- 
tions of a and hence is an increasing function of a. 

It is also worth mentioning that a similar analysis in [9] 
facilitates also to address the issue of successful transmission 
per a given user, and when the average rate per user is 
examined, the results are in agreement with Prop. |6] (or Corr. 

El. 

C. Numerical Results 

The achievable rates of the MCP, SCP and ICFS schemes 
with = 0.5 and or P = 14 [dB], are plotted as functions 
of the erasure rate q, in Figure |6] The figure demonstrate the 
superiority of MCP over both SCP and ICFS schemes. It is 
also observed that ICFS is beneficial over SCP for erasure 
rates smaller than a certain threshold. 

D. Discussion 

A different scenario which gives rise to a similar, but non- 
ergodic, erasure model in a cellular scenario is one in which 



the MTs may be active or not in a given block with a given 
probability accounting for dynamic user activity. According 
to this model, each user is randomly selected to be active 
with probability (1 — g) throughout the entire transmission 
block. This may be implemented by means of control, or 
simply when the MT has no information to send. In these 
cases, the erasures are still i.i.d. along the cells' indices but 
are non-ergodic along the time index for each user With 
optimal MCP, the throughput per active user is an r.v., given 
in the large system limit by 

C'nicp — ua — lini 

^-^^ logdet (/at+i + PHn+iEnE^j^H^j^^^'^ , (25) 

where diag(£Jjv) — (ei, . . . , e„) is a realization of the i.i.d. 
user-activity pattern. Since ^^"^ ^^°°> (1 — g) due to the 
Strong law of large numbers, we can replace tr(£?Ar) with 
{1 — q)N in (l25T l. Following similar argumentations as those 
made in [15], it can be shown that Cmcp-ua converges a.s. to 
(1 — q)~^Cnicp- Turning to the SCP and ICFS schemes it is 
easily verified that the throughputs per active user converge 
a.s. to the corresponding rates of the shadowing setup times 

(i-q)-'- 

In Figure |7] the throughputs per active user of the MCP, 
SCP, and ICFS schemes are plotted as functions of the non- 
active probability q, for — 0.5 and P = 14 [dB]. The 
superiority of MCP over the other schemes is apparent. 
Also we can observe that both MCP and ICFS are robust 
and maintain their throughput constant (for ICFS) or almost 
constant (for MCP) regardless of q, while the throughput 
of SCP drops sharply with decreasing q. This is easily 
explained since neither MCP nor ICFS suffers from inter- 
cell interference, the first due to the fact that it leverages 
all received signals as useful, while the second since it 
avoids interference by frequency planning. On the other 
hand, SCP treats the other cell's signals as noise and hence. 



with more active users, more cells have to decrease their rate 
to overcome the emerging interference. 

VI. Concluding Remarks 

In this work we have considered the input-erasure Gaus- 
sian channel. We have shown that with i.i.d. inputs, the input- 
and output-erasure [5] channels are equivalent. Focusing on 
the two-tap case, we have presented an alternative analysis to 
the one reported in [5]. In particular, in contrast to [5] where 
the rate expression involves a fixed-point equation, the new 
expression is given in the form of an infinite sum, which en- 
ables simple calculation of arbitrarily accurate lower bounds. 
Moreover, the analysis is based on simple properties of tri- 
diagonal matrices and also holds for erasures with memory. 
An alternative simpler high-SNR rate power offset expression 
in the form of an infinite sum has also been presented, 
which contrasts to the one reported in [5] which involves 
a fixed-point equation. Numerical results have demonstrated 
the deleterious effects of increasing erasure rate and reveal 
the nature of the worst-case filter (coefficients with equal 
amplitudes) and that of the best-case filter (either coefficients 
is set to zero). The results have then been used to demonstrate 
the benefits of MCP over SCP schemes for a simple cellular 
uplink with shadowing fading or with user-activity model. 
Other settings that include independent shadowing to each 
BS, and more complex user-activity models, are currently 
being studied. 



Appendix 



We need to calculate 



lim 

N^oo 



E 1 logdet (Jat+i + PHn+iEnE^j^H^^^^ 

Using (|5]l and defining the N x N matrix 

Gn — E^nH^n+iHn+iEj^ , 
we can rewrite ( |26] l as 

C= lim E ( 4logdet(Ijv + PGAr; 

N^oo \ N 



(26) 



(27) 



(28) 



Straightforward calculations reveal that Gat is given ex- 
plicitly by 

en(|/io|^ + l^il^) n = m 

enSn+ihohl n = TO - 1 ,^„, 

N\n,ni = S i , -, ' (2y) 

e„_ie„/i^/ii n = TO+l 

otherwise 



\G 



e.g., for the special case of = 3 we have 

/ ei(|/ioP + eie2hohl \ 

eiea/ij/ii ea (|/io P + I'll |^) e^e^hohl \ . 

(30) 

It is observed that Gat is a tri-diagonal matrix (due to the 
two-tap filter), and that an erasure of the nth input symbol. 



i.e., e„ = 0, sets the nth row and the nth column of Gat to 
zero. Hence, each erasure splits Gat into two block-diagonal 
sub-matrices. The latter is easily explained by the following 
example. Assume = 10 and an erasure realization e' = 
(1, 1, 0, 0, 1, 1, 1, 0, 1, 1). In this case we have 

-f 10 + PGio = 
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(31) 



where a is defined in d?), and c = P/iJ/ii. Examining (ISTT l 
it is evident that 

logdet (Jio + PG10) = 

2 log det {1 2 + PG2 ) + log det ( J3 + PG3 ) , (32) 

where G„ = hI^_^_^H„+i is a tri-diagonal n x n Toeplitz 
matrix (equal to G„ with no erasures). Turning to the general 
case 

logdet {In +PGn) = 

N 

^n(e') logdet (7„ + PG„), (33) 

n=l 

where n(e') is the number of n-length consecutive non- 
erasure sub-sequences in e'. Using (|33] |. we can express the 
rate for any block size N as 



i?/v =E 



N 



logdet (7„ + PG„) 



N 



n(e) 



(34) 



logdet (/„ + PG„ 



Defining l^g. ^.m] an indicator function of the event that 
an isolated sequence of n non-erasures starts at index rn in 
e we get that 



n(e) 



N 



= E 



1 ^ 



m} 



= 1- 



(35) 

where the last equality is due to the fact that the erasures are 
i.i.d, henc^ 



IE [l{e: n,m}] — 



m = 2,...,N -n 







otherwise 



(36) 



-For simplicity matere we neglect the fact that for m = 1, 
E [l{e: n,l}] = 9(1 ~ q)" ™d claime it has no effect in the large system 
Umit. 



Next, we calculate det (J„ + PGn) ; n = 1, . . . , A^, by 
recalling the fact that the determinant of any tri-diagonal 
n X n matrix V obeys the following recursive relation 

det Vn = [Vn]n,n det Vn-1 

det Vn-2 ; n>3 , (37) 

where Vn-i, Vn-2 are the first two main principal sub- 
matrices of Vn (composed of the first n — 1 and n — 2 rows 
and columns of Vn, respectively). Applying the recursive 
formula to £)„ = J„ + PGn we get the following difference 
equation: 

det (Dn) = a det (-D„_i) - 6^ det {Dn-2) ■ n>2 , (38) 

with initial conditions dot (£)o) = 1 and dct(£)i) = a, 
where a and b are defined in Q. A solution to a difference 
equation of this form is given by 



det (£>„) = cpr" - c^s"" , 



(39) 



where r and s are defined in (O. The initial condition are 
used to calculate and 0: 



(p ~ and (f) = 



(40) 



hence, 



det (£>„) 



■ n = l,...,N . (41) 



Substituting (l35]l and (STJ into (l34]i we get that 



N 



N 

(42) 

It remains to show that Rjsi — ~^°°> R. To do that we 
notice that 

1 ^ 

i?iv = - ^ ^ (1 - <7)"n log det - 



n=l 



n=N+l 



Hence, i?jv can be bounded by 



N 



n=l 



n=N 



(44) 

This is since the summands in (l43T l are positive, and by 
noticing that log det (-Dn) < n(i (for some positive /? > 0). 
After some algebra we get that 

R - - /3(1 - qf {{N + l)g + 1)) < i?jv < . 



Nq 

Finally, taking N ^ 00 completes the proof. 



(45) 
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